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ON THE APPLICATION OF TRANSONIC SIMILARITY RULES 
TO WINGS OF FINITE SPAN! 


By Jonn R..SPrprrer 


SUMMARY 


The transonic aerodynamic characteristics of wings of finite 
span are discussed from the point of view of a unified small 
disturbance theory for subsonic, transonic, and supersonic 
flows about thin wings. Critical examination is made of the 
merits of the various statements of the equations for transonic 
flow that have been proposed in the recent literature. Tt is 
found that one of the less widely used of these possesses con- 


siderable advantages, not only from the point of view of a priori ` 


theoretical considerations but also of actual comparison of 
theoretical and experimental results. The similarity rules and 
known solutions of transonic-flow theory are reviewed, and 
the asymptotic behavior of the lift, drag, and pitching-moment 
characteristics of wings of large and small aspect ratio is 
discussed. It is shown that certain methods of data presentation 
are superior for the effective display of these characteristics. 


INTRODUCTION 


The small perturbation potential theory of transonic flow 
proposed apparently independently by Oswatitsch and 
Wieghardt, Busemann and Guderley, von Kármán (refs. 1 
through 6), and others is now supplying a fund of information 
regarding transonic flow about aerodynamic shapes. Solu- 
tions have been given for two-dimensional flow around 
airfoils at both subsonic and supersonic speeds in papers 
by Guderley and Yoshihara, Vincenti and Wagoner, Cole, 
Trilling, Oswatitsch, Gullstrand (refs. 7 through 16), and 
others. In the application of these results to specific ex- 
amples, two items of theoretical interest have been noted 
(see, in particular, refs. 8, 17, and 18): (a) The theoretical 
results appear to be applicable at Mach numbers far removed 
from 1 even though, in most cases, the results have been 
obtained from equations valid only in the immediate neigh- 
borhood of sonic speed. (b) In the application of theoretical 
results to specific examples at Mach numbers other than 1, 
it has been noted that certain ambiguities exist in the 
theoretical determination of aerodynamic quantities. It is 
one of the purposes of this report to investigate these two 
points. This is accomplished by examining transonic flow 
from the point of view of equations that are valid throughout 
the Mach number range rather than only in the neighborhood 
of sonic speed. Such an approach emphasizes the relation 
between the roles of linear theory and of nonlinear theory 
in the transonic range. 


The similarity rules provided by the theory (refs. 5, 6, 


and 19 through 22) have also proved to be useful in the cor- 


relation and interpretation of experimental data. It is 
with the latter aspect of the transonic-flow problem that 
the present paper is primarily concerned. In this paper, 
the similarity rules and their application to the specific 
problem of concise presentation of lift, drag, and pitching- 
moment characteristics of wings are given in detail. The 
known solutions of two-dimensional transonic flow are re- 
viewed and the asymptotic behavior of the aerodynamic 
characteristics of wings of large and small aspect ratios is 
examined. It is shown that certain methods of data 
presentation are advantageous for displaying these charac- 
teristics. 


SYMBOLS 
A aspect ratio ' 
A: [U klti] 2A 
a speed of sound 
Qo speed of sound in the free stream 
a* critical speed of sound 
b wing semispan 
Co drag coefficient 
a (Uk)! 
i Gaye © 
Co, drag coefficient of symmetrical nonlifting wings 
Fi (Wok) 
Do (tje) “Pe 
ACp contribution to drag coefficient due to lift 
(ACp/a*) — [Uk(t/e)** (ACp/ a2) 
CL lift coefficient 
(Czla) [Uke (t/e)}** (Cia) 
On pitching-moment coefficient k 
(Cm/a) [Uk (t/e)]'* (Cala) 
C; pressure coefficient s 
a 1/3 
Cp ee Cs 
Co critical pressure coefficient 
lF. _center-of-pressure function 
c wing chord 
Ca, section drag coefficient of symmetrical nonlifting 
` ` airfoils 
~ (Usk)? ys 
Cag, (H “40 
Ace contribution to section drag coefficient due to lift 


(Aca/ a ) 


[Tke(tle)}* (Acao?) 


! Supersedes NACA TN 2726 entitled “On the Application of Transonio Similarity Rules,” by John R. 8preiter, 1952. ~ 
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ĉi section lift coefficient 

(eda) [Uk te) (ea) 

D drag function i 

D, drag function for symmetrical nonlifting wings 

Da drag due to lift function 

do section drag function for symmetrical nonlifting 
airfoils 

ds section drag due to lift function 

k coefficient of nonlinear term of differential 


equation for yg. (See eqs. (7), (18), (19), (20), 
(23), and (35).) 
lift function 
section lift function 
pitching-moment function 
section pitching-moment function 
local Mach number ` 
free-stream Mach number 
pressure function . ` 
stretching factors defined in equation (B8) 
maximum thickness of wing 
free-stream velocity 
V velocity components parallel and perpendicular, 
i respectively, to the flow direction -ahead of a 
shock 
Cartesian coordinates where x extends in the 
direction of the free-stream velocity 
distance from wing leading edge to center of 


pressure 


a 


AST. RESSA 


angle of attack 


a 
(t/c) 
y ratio of specific heats, for air y=1.4 
A arbitrary constant 

(M—1) 
: [Tok Ce)” 
(M— 1) 
Si VAO 
T ordinate-amplitude parameter 
i) velocity potential 
p perturbation velocity potential 


SUBSCRIPTS 


l values given by linear theory 

W conditions at the wing surface _ 

J conditions immediately upstream from shock 
conditions immediately downstream from shock 


FUNDAMENTAL CONCEPTS 


BASIC EQUATIONS 


The quasi-linear partial differential equation satisfied by 
the velocity potential # of steady isentropic flow of a perfect 
inviscid gas can be expressed in the form (ref. 23, p. 25) 


(a?@—,7) Prt (a? — 8,7) ,, + (V— 8,7) ,,— 29,8, 8,4— 


20,9, 0,,—2P,P,P,,—= 0 (1) 


where the subscript notation is used to indicate differentia- 
tion anda is the local speed of sound given by the relation 


ee — -——_— i ee te. 


ordinates of wing profiles in fractions of chord - 
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=at I @2+8)+2— Uo’) (2) 
In this latter equation U, and a, are, respectively, the ve- 
locity and speed of sound in the free stream and y is the ratio 
of specific hests (for air y=1.4). 

Introducing the perturbation velocity potential , where 


y=—U ate (3) 


it is possible to express equation (1) in terms of the deriva- 


tives of ¢ as follows: 


(1 =; o Pex Sy PP: 


e| (y-+1)U e+ E 


pr + 14 lor + oi) | + 


er) (y—1)U set ettet oF |+ 


Pet =D Uet 2S (p+ pE pë |+ 
(4) 
2 a Pol Ust gp) + 


2 a vel U+ Pz) + 
; 2 A pyp: E 


If it is assumed that all perturbation velocities and perturba- 
tion velocity gradients (represented by first and second 
derivatives, respectively, of y) are small and that only the 
first-order terms in small quantities need be retained, equa- 
tion (4) simplifies to the well-known Prandtl-Glauert equa- 
tion of linear theory 


s) Prz F Pyy t Pru ™=0 (5) 


where the free-stream velocity is directed along the positive 
z axis as shown in figure 1 and where M, is the Mach number 
of the free stream. It is well known that equation (5) leads 


(l= 


Z 


FiavurE 1.— View of wing and coordinate system. 


* 
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to useful results in the study of subsonic and supersonic 
flows about thin wings and slender bodies but that it is 
incapable, in general, ot treating transonic flows. The 
failure of linear theory in the transonic range 1s evidenced 
by the calculated value of p, growing to such magnitude 
that it can no longer be regarded as a small quantity when 
compared with U,. 

Second-order theory for thin wings would involve solution 
of the equation 


(1 —M yst Py t Prr: = 
M, ve gati — Phen t Pu) tT U, = (even oem) | (6) 


Actually, we are interested in retaining higher-order terms | 


only to the extent that is necessary to allow the study of 
transonic flow. Examination of the known characteristics 
of transonic flow fields indicates that the first term on the 
right can often become of importance and should be retained. 
The remainder of the terms on the right can never become 
large for transonic flows about thin wings at small angles of 
attack and can be safely disregarded. The simplified equa- 
tion is 


(1— 


0) east Ot r= Me, oa P29 korpa (7) 

It follows from the basic assumptions of small-disturbance 
transonic theory that equation (7) is valid everywhere in the 
flow field, both fore and aft of any shock waves, but cannot 
provide any information on the discontinuities in that 
represent the shock waves. This information must be ob- 
tained through consideration of the classical equations for 


the velocities on either side of an oblique shock. Thus- 


if- U, represents the velocity immediately in front of the 
shock wave and U2 and V2 represent the velocity components 
parallel and perpendicular, respectively, to U; that occur 
immediately behind the shock, the equation for the shock 
polar (ref. 23, p. 108) provides that 


U, U2 ER 


V= ( U,— = 
U, aF U,0,+a* (8) 


Ur)’ — 
yrl 
where a* represents the velocity of sound at a point where 
the local Mach number is unity and is commonly designated 
the critical speed of sound. It is related to the free-stream 
velocity and speed of sound by the following equation de- 
rived from equation (2): 

yet =a’ tI U? (9) 
Except for the important case of the bow wave in supersonic 
flow, U; is not, in general, alined with the direction of the z 
axis, but is inclined a small angle. With the resolution into 
components parallel to the axes of the coordinate system and 
the retention of only the leading terms, in a manner similar 
to that used in the derivation of equation (7), equation (8) 
provides the following relations between the velocity com- 


ponents (potential gradients) immediately fore and aft of 
the shock: - 
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U’ — a*? + U.( + T ] 
OEE E E re S 
y+ ee 


(10) 


- This result can be put into a more desirable form by making 


use of the — relation derived directly from equation 
(9): 
2 Mé,*—-1 


2 es ee 
Z ytd Mè 


U (11) 


Substitution of this expression into equation (10) and rear- 


rangement of the terms yields the desired result 
(1 —M,’) (Pz T Pz) t (m Pr) +(ePz, = Pe) = 


This equation corresponds to the shock-polar curve for 
weak shock waves inclined at any angle between that of 
normal shock waves and that of the Mach lines. 

In addition to satisfying equations (7) and (12), the per- 
turbation potential must provide flows compatible with the 
following physical requirements: (a) The flow must be uni- 
form far ahead of the wind, and (b) the flow must be tan- 
gential to the wing surface. Therefore, the following boun- 
dary conditions are to be oe for the perturbation 
potential: 


at r= — oœ 
(gz) = CARa (p:)o= 0 (1 3) 
at the wing surface W | 
Pz _oZ 
(Tir) r Gy 


where òZ/òzx refers to the local slope in the x direction of the 
wing surface. A systematic application of the perturbation 
analysis indicates that the boundary conditions specified on 
the wing surface should be simplified by approximating the 
fraction on the left by ¢,/U, Furthermore, it is consistent 
with the assumption of small disturbances to satisfy this 
boundary condition on the two sides of the zy plane rather 
than on the wing surface. Equation (14) is therefore re- 


placed by 
(+ e)s = aa Gt oe) 
where the shape of the wing profile is given by 
~ Zle=rf(x/c, y/b) (16) 


where f(z/c, y/b) represents the ordinate-distribution function 
and r is an ordinate-amplitude parameter. Note that, in 
general, a variation of r represents a simultaneous change of 
the thickness ratio, camber, and angle of attack. In the 
special case of a nonlifting wing having symmetrical sections, 
r is proportional to the thickness ratio; for lifting flat-plate 
wings of vanishing thickness, 7 is proportional to the angle 
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of attack. In order to obtain unique and physically im- 
portant solutions, it is necessary to. assume the Kutta condi- 
tion (that the flow leaves all subsonic trailing edges smoothly). 

Upon solving the above boundary-value problem for the 
potential, one may determine the pressure coefficient by 
means of the formula 


=T, Pz (17) 


It should be noted that the results obtained by using the 
foregoing approximate equations might be expected to tend 
toward those of linear theory as the free-stream Mach number 
Me, departs far from unity. This follows from the fact that 
the product pz% becomes small relative to the linear terms 
under this condition. Solutions of the equations for tran- 
sonic flow found to date have all possessed this property. 


COMPARISON WITH OTHER STATEMENTS OF THE TRANSONIC-FLOW 
EQUATIONS 


y 
As a result of minor variations in fhe perturbation analy- 
sis, recent papers have used at least four different relations 
for k, the coefficient of the nonlinear term in the simplified 
equation for the perturbation velocity potential. As indi- 
cated in the preceding paragraphs, straightforward develop- 
ment of the theory leads to the relation 


—w2rtl 
k=M, Uy, (18) 
This is sometimes simplified (e. g., refs. 22 and 24) to 
_rytrl 
k= T, 9) 


by arguing that M, can be set equal to ynity in this term 
without much loss in accuracy, since the right-hand side of 
equation (7) is merely an approximation to allow the treat- 
ment of transonic flows and rapidly diminishes in magnitude 
as M, departs from unity. In some treatments (e. g., refs. 
16 and 21), equation (1) is divided by a’, and the quotient 
l/a? in each term is expanded in a binominal series. When 
this is done, the coefficient k of the term involving gps is 


k= M2 pes] (20) 


Still another expression for k is used by Oswatitsch in 
references 13 and 14. Two derivations are given, one based 
on mass-flow considerations (ref. 14) and the other (ref. 13) 
-based on simplifying equation (1), under the assumption 
of nearly parallel flow, to 


(1—M") vert ogy +¢u=0 (21) 

and expanding the variable coefficient 1—M in the series 
_ ape 1 are 1— Mo Ov 

1—M=1—- M3 SP (22) 


where M is the local Mach number and a* is the critical 
speed of sound as defined im equation (9). Comparison of 
equations (21) and (22) with equation (7) shows the coeffi- 
cient k in this approximation is 


and (12). 
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2 2 
bate (23) 


It should be noted that the four alternative relations for k are 
identical for M,=1 and all but that given by equation (19) 
are zero for M,=0. 

A similar situation arises in the derivation of the simplified 
equation for the shock polar. Here again the precise form 
of the expression for the coefficient k of equation (12) depends 
on the details of the perturbation analysis. The most 
important point from a practical point of view is that the 
same expression for k is used in both the equation for the 
potential and that for the shock polar, namely equations (7) 
While this point has not always been explicitly 
stated, it is actually a necessary condition for the existence 
of the transonic similarity rules. 

Although each of the above alternative forms of the 
perturbation equations has been used at least once in the 
recent theoretical investigations of transonic flow, the most 
widely used set of equations are those derived under the 
more restrictive assumption that all velocities are small 
perturbations around the critical speed of sound a* rather 
than around the free-stream velocity U,. In the latter 
scheme, the perturbation potential is defined by (see, e. g., 
ref. 6 or 20) 


yp’ =—a*r +Ë (24) 
and the resulting differential equation for ¢’ is 
1 
P'wt ett 0's) 2x (25) 


The approximate relation for the shock polar is 


zt 
G no P AG ‘1 ey) = atti (fits £ CEEA (26) 


The corresponding boundary conditions aro specified as 
follows: 


at r =— © 
(y'x)o=U o—a* s= a ~M’), (o'o (p :)o= 0 
at the wing surface 
E EE (2 y 
CORTA EEDU m b (27) 


where the shape of the wing profile is still given by equation 
(16). The equation for the pressure coefficient is approxi- 
mated similarly, thus, 


0y = Alo lodd (28) 


This statement of the equations for transonic flow is clearly 
identical to those given previously, herein, when the free- 
stream Mach number is unity. Although the derivation of 
the a* equations requires that the free-stream Mach number 
be very close to unity, these equations have been used with 
good success by a number of authors to calculate the aero- 


dynamic forces on airfoils at Mach numbers considerably 
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removed from unity. In so doing, it has been suggested that 
it might be preferable to use more accurate relations for the 
pressure coefficient or the boundary conditions; for instance, 
it has been suggested that a* be replaced with U, in the 
equation for C,’. This matter has been discussed at length 
in references 8, 17, and 18. Since no restriction requiring 
the Mach number to be near unity is made in the U, analysis, 
it is informative to examine the relation between the results 
of the a* and the U, analyses. ~This is done in Appendix A. 
It is found that the a* analysis, if performed in a completely 
consistent manner using equations (24) through (28), yields 
values for C, that are identical to those given by the more 
general U, analysis using k=(y+1)/U,. This somewhat 
paradoxical result is achieved through the action of a number 
of compensating effects and only applies to the pressures and 
the forces and moments derivable therefrom. It should be 
noted, in particular, that the values of the local velocities 
and Mach numbers provided by the a* analysis for flows 
having free-stream Mach numbers other than unity are in 
error. Throughout the remainder of this report, the discus- 
sion will be based on the U, analysis. 

A significant case where the alternative relations for k lead 
to different results is the prediction of the critical pressure 
coefficient Cp., defined as the value of the pressure coefficient 
O, at a point where the local Mach number is unity. It is 
important that a reasonably good approximation be main- 
tained for the variation of Cp„ with M, because shock waves 
make their first appearance, and the airfoil first experiences 
a pressure drag when C, becomes more negative than Cp, 
somewhere on the airfoil surface. In the present approxima- 
tion, Cp„ corresponds to that value of O,, and, hence, or 
Pz, at which equation (7) changes locally from elliptic to 
hyperbolic type. This condition is recognized by the vanish- 
ing of the coefficient of pzs, thus 


1—M—k(¢2) --=0 


or, in view of equation (17) 
t 


2 2 
O=- Fr (e=— Er -MA (29) 


The exact relation for isentropic flow is (e. g., ref. 25, p. 28) 


2 |p 2 yl a Vr 7 
Ore PERE +) 1 ow 


The variation of C,_ with M, has been computed by use of’ 
the exact relation and each of the four approximate relations. 
The results are presented graphically in figure 2. It may be 
seon that a reasonably good approximation for C,_ is ob- 


tained over a wide Mach number range when k is taken as~ 


given in either equation (18) or (23), and that a somewhat: 
greater error is incurred when equation (20) is used. On the 
other hand, equation (19) leads to a very poor approximation 
for Cp... 


Similar comparisons can be made for local Mach numbers 
M other than unity by noting that the coefficient 1—142—ke, 
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: -8 
-Me [2+(y-1)ME]/Up 
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Fieurb 2.—Variation of critical pressure coefficient with Mach number. 


of yr; m equation (7) corresponds, in the present approxima- 
tion, to 1— AP, thus 


k 


1—M=1—M?—kp=1—Mè +0, (81) 
The corresponding exact relation for isentropic flow is 
P 1+4 Me pa 
o=] -1+ — (32) 


147-1 ap 
2 i 
The results so obtained are generally similar to those in- 
dicated in figure 2, although the relative accuracy of the bet- 
ter approximations changes somewhat with the situation. 
All the approximations are exact, of course, when C,=0; on 
the other hand, none of the approximations are exact, except 
for isolated cases, when C, is different from zero, even though 
all the approximations agree among themselves when the 
free-stream Mach numberis unity. In order to provide some 
information regarding the errors that are likely to be in- 
curred when C, is not very small, figure 3 has been prepared 
illustrating the variation of local Mach number with pressure 
coefficient for a free-stream Mach number of unity. 
A second case where the exact and approximate relations 
may be compared is furnished by considering the velocity 


- =— -— = See fe pe 


D 2 


~ 


Js 
Approximate 


1.0 > 


Fiaurn 3.—Variation of local Mach number with pressure coefficient 
for M,=1. 


jump through a shock wave. If the flow abead of the shock 
wave is uniform and parallel to the x axis, the results may be 
conveniently represented by the shock-polar diagram in 
which Vz/a* is plotted as a function of U2/a*. The exact 
relation is furnished by equation (8). The corresponding 
approximate relations are determined from equation (12) by 
setting Px) Pry» Py, 2nd pz, equal to zero, whereby U,=U,, 
M,=M,, and 


k 

e.2-|— MÀ t5 Pn | Pary (33) 
Once the variation of ¢., with p; is determined for a given 
Me the corresponding variation of Va with Ue may be 
readily determined since, for this case, 

Us=Uo+ Pays (34) 
Va 
a* 
using both the exact and approximate relations, and the 
results are presented graphically in conventional shock- 
polar form in figure 4. This figure contains, in addition to 
curves for the four expressions for k discussed previously, a 
curve computed using 


= y+1 
k=M,.* T. 


V= Pre 


The variation of — with 5 for M4,=1.2 has been computed 


(35) 


i/o 


-Mo [2+(y-Me | /U, 


jf Mo (7 t0U 


Figure +—Exact and approximate shock polars. 


- 
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This expression for k arises in a series expansion of the right- 
hand member of equation (10). (For sake of completeness, 
the corresponding variation of C,,. with Me is also included 
in figure 2.) Just as with the comparison of the critical 
pressure coefficient, the expressions for k given by equations 
(19) and (85) lead to poor approximations, whereas any other 
of the expressions lead to reasonably good approximations. 
A notable point is that the expression for k given by equation 
(18), that is k=Mè(y+1)/U., leads to the exact relation for 
velocity jump through a normal shock wave. 

In order to facilitate comparison with previous results and 
to achieve an economy of notation, the present analysis is 
carried as far as possible without specifying a particular rela- 
tion fork. That is, the equations of the analysis and reduced 


` parameters with which the results are expressed are written 


containing k which may be equated to any of the five stated 
relations, or in fact, to any coefficient that does not depend 
on x,y,zor . The actual values of the pressure coefficient 
and Mach number for an airfoil of specific thickness ratio, 
however, depend on which relation is selected fork. Whore- 
ever such values are given, they will be those obtained by use ° 
of the expression for k given in equation (18), that is, 


are Ytl 
k= M, Uy, (18) 
The principal reason for this choice is that it appears to 
provide a set of equations, or a mathematical model, which 
approximates certain essential features of transonic flow with 


` superior accuracy. 


RELATION BETWEEN TRANSONIC THEORY AND LINEAR THEORY 


It is important to recognize that wing theory based on 
equation (7) is valid for all Mach numbers below the hyper- 
sonic range. At subsonic and supersonic speeds, equation (7) 
is of the same order of accuracy as the Prandtl-Glauert 
equation of liner theory (eq. (5)) although more difficult to 
solve. At -M.=1, equation (7) is identical with equation 
(25), now widely used in the study of transonic-flow problems. 

On the other hand, there is no a priori method for de- 
termining whether or not a solution of the equations of 
linear theory will be valid in the transonic range. One can 
only decide by solving the problem under the assumptions of 
linear theory and then inspecting the magnitudes of the 
terms, particularly of ¢,, to see whether or not they can be 
regarded as small quantities. If the terms are sufficiently 
small, the linear-theory solution is presumed valid oven 
though the Mach number may be near unity. Linearized- 
theory solutions have been obtained for a great number of 
practical wing problems and their behavior in the transonic 
range is now well known. To review briefly: For unswept 
wings of infinite span, linear theory indicates that the 
magnitude of s on the surface of a given airfoil is propor- 
tional to 1/-/[I—14,3|; consequently, z approaches infinity 
as M, approaches unity and the theory is clearly inapplicable. 
For wings of finite span, however, the perturbation velocities 
may be large or small at sonic velocity, depending on tho 
particular problem as discussed in detail in reference 26. 
Specifically, for three-dimensional lifting surfaces of zero 
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thickness, the velocities remain finite everywhere except at 
the leading edges, their magnitudes generally increasing with 
increasing aspect ratio and angle of attack. For wings of 
nonzero thickness, however, s generally becomes large 
logarithmically as 1—M,? approaches zero;. consequently, 
linear theory is inapplicable within some Mach number 
range surrounding unity. 

Summarizing, linear theory is applicable to lifting surfaces 
of small or moderate aspect ratio at all transonic speeds, 
but fails for wings of finite thickness within a range of 
Mach numbers surrounding unity. The range of inapplica- 
bility diminishes to zero as the aspect ratio, thickness ratio, 
and angle of attack of the wing tend to zero. 

In treating transonic flows for which linear theory is 
applicable, it is often advantageous to consider the special 
case of sonic flow (M,=1) separately. Equation (5) for 
the perturbation potential then reduces to a particularly 
simple form i 

Pyt Pz = 0 


Solutions of this equation, in conjunction with the boundary 
conditions given by equations (13) and (15), are identical to 
those of linear theory found by solving equation (5) and 
subsequently setting M,=1, but can be obtained with much 
less effort. Since, in addition, the results of this simple 
theory, now generally known as slender-wing theory, are 
also applicable to low-aspect-ratio lifting surfaces throughout 
the entire Mach number range, a considerable number of 
solutions of slender-wing theory have been presented in the 
last few years (e. g., réfs. 27, 28, and 29). These results 
are, of course, applicable to flows at M,=1 to exactly the 
same extent as the results of linear theory. 


SIMILARITY RULES 


In reference 6, von Kármán derived similarity rules for 
the pressure distribution, lift, drag, and pitching moment of 
airfoils in transonic flow using equations (24) through (28). 
The same equations were used in reference 20 to determine 
the transonic similarity rules for wings of finite span. The 
corresponding similarity rules of linearized subsonic and 
supersonic wiùg theory were also derived’ and compared 
with the transonic similarity rules in the latter reference. 
It was shown that the similarity rules of linear theory con- 
tain an arbitrary parameter and can be expressed in many 
forms, one of which coincides with the similarity rules of 
transonic flow. 

A derivation of the transonic similarity rules, based on 
the U, equations with unspecified k, is provided in Appendix 
B. This derivation possesses the advantage of being based 
on a single statement of the problem of wing theory that 
is uniformly valid at subsonic, transonic, and supersonic 
speeds. It follows from the preceding discussion that the 
results so found are identical to those of references 6 and 
20 if k is equated to (y+1)/U,. The similarity rules for 

Op, Cr, Cm, and Cp are given in Appendix B as follows: 
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om giin ?[Gie (Tak? VIL? A; > z| (36) 


2/3 1—M, 
Oaa <| oye TE ra] k 
2/3 1 —M,* 
Comer EA] es 
5/3 L=M; 
Co= EB D (ks J- 7s vl —M, 24 | (39) 


where the geometry of related wings is given by equation (16): 
a (Z/e)=rf(e/c, y/b) 


Equations (36) through (39) are functional equations. For 
example, equation (36) is to be interpreted as stating that 
the pressure coefficient C, is equal to 774/(U,k)'* times some 
function of a number of specified parameters. The 
foregoing equations have been Uea for flows where M, <1. 
If M,>1, the radical Y1—M,? should be replaced with 

Mê—1. The functions P, G M, and D are different, 
however, for subsonic and supersonic flow. Consequently, 
subsonic flows may be related to other subsonic flows by 
the similarity rules, but not to supersonic flows, and con- 
versely. 


ALTERNATIVE FORMS OF THE SIMILARITY RULES 


It is important to recognize that the similarity parameters 
may be combined or regrouped in any manner whatsoever, 
provided the same number of independent parameters is 
always retained. For instance, in much of what follows, it 
will be found desirable to use the square of -~1—M2/(Ukr)'* 
and to replace ¥1—M,7A with a new parameter (Ukr) BA: 
obtained by dividing y1—M ŻA by ¥1—M2/(U,kr)!*. In 
terms of these parameters, ihe similarity rules are 


o= Guan P UE (Ukr)? A; 2, zg] (40) 
=op |a Ukya] an 
C= M aye (Tok A] aa 
Co= cra D| aa (Wyler) A] (43) 


wherein the geometry of related wings is again given by 
equation (16). 

The similarity rules thus formulated are totally equivalent 
to those given by equations (36) through (39) but possess 
three outstanding advantages: 

(a) The indeterminacy at M,=1 resulting from two 


| parameters simultaneously vanishing is removed. 
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(b) The squaring of the first parameter avoids the neces- 
sity of changing parameters as sonic speed is passed. 

(c) The use of the parameter (U,tr)’* A rather than 
J1-M A aids in distinguishing the regimes in which linear 
theory is applicable in the transonic range from those in 
which nonlinear theory must be used. Thus as (U,kr)* 
A approaches zero, linear theory is always applicable, pro- 
vided, in some cases, that M, is not precisely equal to unity. 
On the other hand, as (Ukr) A becomes large, linear 
theory is not applicable in the transonic range and nonlinear 
theory must be used. 

The forms of the similarity rules given in equations (36) 
through (43) have been the source of some confusion, due 
to the multiple role that 7 plays in determining the thickness 
ratio, camber, and angle of attack. As can be seen from 
equation (16), all three of these geometrical quantities are 
linearly proportional to r. A more explicit statement of the 
symmetrical profiles of nonzero thickness may be obtained 
by rewriting the expression for Z/c in terms of the thickness 
ration t/c and angle of attack a rather than 7, thus, 


Zle=tle| 9 E “e z] 


Comparison of equations (16) and (44) indicates that t/c 
plays a similar role to r but that it is necessary to introduce 
a second parameter aj(t/c). In this way, a new set of equa- 
tions expressing the similarity rules is obtained which 
correspond to equations (40) through (43), although expressed 


(44) 


in terms of a and t/c rather than 7. They are 
6, Fede O= Plinii; sle, 418) (45) 
GC Co Lb A, 2) (46) 
OnT 0 MG, 4,3 (47) 
n=l Co=Dlb, 4,8) (48) 


where £, A, and @ are similarity parameters defined ag 
follows: 


_ M-il 
E= A 


SLOPE OF PRESSURE CURVE AT Mo=1 ot 


Ä=KUk) tA, =% 


Ho (49) 


Liepmann and Bryson (refs. 17 and 18) have made the 
following observations which enable the determination by 
simple and intuitive considerations of the slope of the ©, 
versus MM, curve at M,=1. It is a well-known fact that, at 
slightly supersopic Mach numbers, the detached bow wave 
is far away from the airfoil and nearly normal. It is also 
well known that the Mach number downstream of weak 
normal shock is as much below unity as the Mach number 
upstream is above unity. Consequently, the Mach number 
distribution on the airfoil should be independent of Mach 
number in the neighborhood of M,=1, that is, 
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(50) 


(TIE) ss 


Now, if the isentropic flow relation for C, in terms of M and 
M, 


L uM? ri 


2 


C= i (51) 


is differentiated with respect to M., Me is equated to unity, 
and ‘equation (50) is introduced, the following relation, first 
given in reference 8 by Vincenti and Wagoner, is found: 


(am, Mo=1 “Hi “oy (Co) moat (52) 


Since the above-mentioned derivation is based to a certain 
extent on physical reasoning and makes use of the exact 
rather than approximate relation between pressure and Mach 
number, it is of interest from the present point of veiw to 
review the equivalent result contained in the model of tran- 


` sonic flow provided by equations (7), (12), (13), (15), and 


(17). In common with the other characteristics of transonic 
theory, the result can be expressed conveniently in the form 
of a similarity rule. Thus, recall the approximate relation 
for the local Mach number given by “ae (31). 
1—-M?= 


1—M ag (31) 


Now, aie’ to the similarity rules, transonic theory 
does not provide information about C, or M, alone, but only 


about parameters such as C, and >» Thus, the similarity 
rule for local. Mach number is the following: 


M-i _ Mè- Oy _, 
[= O Or ayr CTE O (63) 


-The solutions of the equations for transonic flow obtained 


for slightly supersonic flow by Vincenti and Wagoner in 
reference 8 and for slightly subsonic flow by Cole in reference 
11 indicate that the approximate relation which corresponds 
to the Mach number freeze is the following: 


GE | (64) 


It may be recognized that this relation is equivalent to the 
exact relation given in equation (50) if k is independent of 
M., as im equation (19). This is not the case, however, 
when the presently preferred relation for k, namely, equation 
(18), is used. Given equation (54), differentiation of equa- 
tion (53) shows that the slope of the pressure curve in the 
reduced parameters is 
dc, 
2) = 2 (55) 


dg, 


It is a simple matter to derive the corresponding slope for 
the C, versus M, curve at M,=1, provided a specific expres- 
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sion is selected for k. If kis equated to M,7(vy+1)/U., the 
final result is 


2 
(ae. M,=1 =i 3 (Co) samt 656) 
It is interesting to note that if k is independent of Ma, the 
slope of the pressure curve is 


(a) rH 67 


Thus, although equation (19) leads to the exact relation for 
the rate of change of the local Mach number, the slope of 
the pressure curve is considerably less accurate than that 
provided by equation (18). 

The range of problems tó which the foregoing results 
apply is not known at present. Inasmuch as the entire 
phenomenon. appears to be connected with the presence of a 
detached bow wave which stands normal to the flow, intuitive 
considerations suggest that the results are probably appli- 
cable at least to symmetrical airfoils at zero or infinitesimal 
angles of attack but perhaps not to airfoils at larger angles 
of attack or to wings of finite span. 


_ 


APPLICATIONS 
FUNDAMENTAL HYPOTHESES AND PRINCIPLES 


The remainder of this report is principally concerned with 
the deduction of the qualitative, and to some extent quantita- 
tive, characteristics of thin wings in transonic flow by.«means 
of simple logical considerations based primarily on the simi- 
larity rules together with the following hypotheses: 

(a) Nonlinear theory based on equation (7) is applicable 
to all problems. 

(b) Linear theory based on equation (5) is valid for all wings 
at Mach numbers either appreciably below or above unity. 

(c) Linear theory is valid at all Mach numbers, except 
possibly very near unity, for wings of small aspect ratio. 

(d) The differential pressures between the upper and lower 
surfaces of a wing having symmetrical airfoil sections are 
proportional to the angle of attack for at least a small range 
of angles about zero. 


(e) The slope of C, versus $, at M4,=—1, defined by equation 
(55), is applicable at least to symmetrical airfoil sections at 
zero or infinitesimal angles of attack. 

The consequences of the foregoing statements will be con- 
sistently pursued in the following sections in the discussion 


of the aerodynamic characteristics of airfoils and complete 


wings. Throughout, the analysis will be restricted to wings 
having symmetrical profiles. Whenever specific results are 
to be used to illustrate the statements, they will nearly always 
be for symmetrical-wedge or double-wedge profiles and for 
wings of triangular plan form. This choice is dictated by the 
availability of theoretical results. 

The basic principle in the following analysis is to express the 
similarity rules in such forms that the lift, pitching-moment, 
and drag coefficients can be studied for limiting values of the 
parameters with no chance for ambiguity due to indetermi- 
nate forms. In this respect, the statement of the similarity 
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rules provided by equations (45) through (49) will be found 
particularly useful. 


PRESSURE DRAG OF SYMMETRICAL NONLIFTING WINGS 
The similarity rule for the pressure drag coefficient of sym- 
metrical nonlifting wings having profiles given by 
(Z/e=(/c) g(ælc, y/b) (58) 


is obtained from equations (44) and (48) by setting a, the 
angle of attack, to zero 


T CU ok)¥8 


C= (oy? Co, =D. (oå) om) 
where 
E= Ee A= [Wok WON A (49) 


‘Therefore, drag results for symmetrical nonlifting wings 


should be presented by plotting the variation of the general- 


ized drag coefficient TA with £, and A. 


At Mach numbers sufficiently removed from unity for lin- 
ear theory to apply, Cp, must be independent of k since + 
does not appear in either the differential equation or bound- 
ary conditions of linear theory. This implies that the param- 
eters £, and A must be arranged in such a manner that & is 
canceled completely from the above relation for Cp,. The 
only alternative inside the function D, is to form the product 
Jt, 4. Numerous possibilities exist on the outside of D, 
depending on whether k is canceled by using ¢,, 4, or some 
combination thereof. Although any of these are acceptable, 
the first will be preferred. In this way the following results 
are obtained: 


Z (j) ¢ |M?—1| _ Aig \M2—1| y 
(CDe): (Uk)? UE Ta} Dal rece (tJe)? ? 
(7) Yop? A |= “TET Da UMETA 


where it should be recalled that D,, is a different function for 
subsonic and supersonic flow. Equation (60) is equivalent 
to the extended Prandtl-Glauert rule. For subsonic flow, 
D’Almbert’s paradox edn that the drag be zero; there- 
fore, for all wings, 


(Co): =0, 


(60) 


—— 
(Co,),=0 (61) 
E.<o 

For supersonic flow, wave drag exists which depends on 
¥M,—1 A as well as on the plan form and airfoil section. 
The general functional relation for the drag coefficient of a 


family-of affinely related wings at zero angle of attack, as 
given by linear theory, is 


(C2), -DNM M#iA), (On, h be Do (hi? A) 


(62) 


Wings of infinite aspect ratio.— For wings of infinite span 
(or airfoils), equation (62), representing the functional 


te eu Sr RNS yee PN ye, ce ggg M e a eee ee 


a e a M > e 
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relation of linear theory for the drag coefficient, reduces to 
the following: 


(ca) = TRX const., (C4 ie comse (63) 
M.>1 


where the value of the constant depends on the shape of the 
airfoil. Numerous experimental data show variations con- 
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Since calculations have been made of the drag in transovic 
flow of simple symmetrical sections at zero angle of attack, 
it is not necessary to speculate further regarding the vari- 


ation of a witb >» At present, complete theoretical infor- 


mation exists for the drag of symmetrical double-wedge 


airfoils throughout the transonic range. Solutions for this 
section have been obtained by transforming the nonlinear 


sistent with equation (63) at Mach numbers greater than | differential equation of the small disturbance transonic 


about that of shock attachment. At Mach numbers closer 


theory into hodograph variables and taking advantage of 


to unity, however, the theoretical values provided by this | simplification of the normally difficult problems relating to 


equation are unreliable. It is evident from inspection of the 
results of linear theory itself that such a failure occurs, since 
the perturbation velocities assumed to be small in the deri- 


the boundary conditions by restricting attention to polyg- 
~ onal profiles. The results for flows having subsonic, sonic, 
and supersonic free-stream velocities have been given, 


vation of the equations are found to become infinitely large | respectively, by Trilling (ref. 12), Guderley and Yoshihara 


as the Mach number approaches unity. . It is apparent, 
therefore, that it is mecessary to resort to nonlinear theory 
for the calculation of the drag of airfoils in the transonic 


speed range. 


A similarity rule for the generalized section drag coefficient - | 


of symmetrical nonlifting airfoils which is valid throughout 
the Mach number range may be obtained from equation (59) 
by setting A equal to infinity. 


Ue | _ Tok" 
Toy > oy 


At a Mach number of unity, the similarity parameter £, 
vanishes and the function ¢,(&,) is a constant 


(Co,) = Dib 2)= =d LE) (64) 


(ca) =d,(0)=const., (caa) = (Uk ian X const. (65) 


indicating that the section drag coefficients of affinely related 
airfoils are proportional to the five-thirds power of their 
thickness ratios. If hypothesis (e) is accepted, the variation 


of ti. with & at A4,=—1 is found to be zero for completed 


airfoils 
dea, ac(Z/t)| 3 ad(Z/t) z Tz 
=) =f (te Te) Emo dale) | suf 


If k is taken to be AZ,?(y+1)/U., (dC,/dM],) 1,1 is given by 
equation (56) and the slope of the drag curve is 


dea, 
dM],/a1,=1 =0( P- =) 


=f zH- — 3 ae ol 


The corresponding exact relation can be determined simi- 
larly by use of equation (52) for the slope of the pressure 
curve. Itis 


ae M,=1 A FHT Ou z(a 


O= = 


0 (66) 


2 
—3(¢e,) Ma=1 


(67) 


-H (Cao) nee (68) 


(ref. 7), and Vincenti and Wagoner (ref. 8). The linear- 
theory solution for pure supersonic flows has been given by 
Ackeret (ref. 30). All of these results are combined on a 
single graph in figure 5. It may be seen that the preceding 


Liepmann and Bryson—-—e -—— -Guderley and Yoshihara 


\ - Vincenti and 


Trilling- Wagoner 


——— Nonlinear theory 
~=- Linear theory 


Figure 5.—Theoretical drag of double-wedge airfoil. 


|. remarks concerning the relation of the results of linear 


theory and nonlinear theory, and the slope of the drag curve 
at a Mach number of unity are illustrated by this comparison. 

Osawatitsch has given approximate solutions for M,<1 
for the pressure distribution on symmetrical biconvex 
airfoils (refs. 13 and 14) and for the pressure distribution 


and drag on NACA four-digit symmetrical airfoils (ref. 14). 


This work has recently been extended to several NACA 
6-series airfoils by Gullstrand (refs. 15 and 16). Their 
drag results are generally similar to those indicated for the 


-double-wedge section in figure 5. 


In problems such as we are considering herein, the final 
test 1s provided by comparison with experimental results. 


‘Unfortunately, experimental results for the transonic speed 


range are scarce, as well as difficult to obtain, and no data 
are available. for direct comparison with the theoretical 
results summarized in figure 5. At present, however, com- 
plete information, both theoretical and experimental, does 
exist for a single-wedge section followed by a straight section 
extending far downstream. In accord with some of the ` 
original papers on this subject, the single-wedge section is 
considered as the front half of a symmetrical double-wedge 
airfoil having a chord c. Solutions for this section obtained 
using transonic flow theory have been given for flows having 
subsonic, sonic, and supersonic free-stream velocities, re- 
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spectively, by Cole (ref. 11), Guderley and Yoshihara (ref. 0); 
and Vincenti and Wagoner (ref. 8). These results are shown 
in figure 6 together with the corresponding experimental 


(y+!) 1/3 


Wess f0 


T bag 
i TT ¥ 


[iy +1) (7/c 173 
Fiaure 6.—Drag correlation of single-wedge section, k=(y+1)/U.- 


results of Liepmann and Bryson (refs. 17 and 18). Both 
the theoretical curve and the experimental points are plotted 
in the same manner as presented in the original papers. As 
remarked in a preceding section, the theoretical results are 
based on a set of equations equivalent to the present equa- 
tions with & equated to (y+1)/U,. The small vertical lines 
on the experimental data points represent the uncertainty 
of the values. The slope of the drag curve at M,=1 is no 
longer zero as indicated for the complete airfoil by equation 
(66) but takes on a positive value given originally by Liep- 
mann and Bryson (ref. 17) and readily derivable from 


equation (66). 
te) 
dé Eno 


This figure indicates that the theoretical and experimental 
results are only in general qualitative agreement. Since it 


(69) 


is shown in a preceding section that theoretical considera- ` 


tions suggest that a superior theory results if k is equated 
to M,/7(v+1)/U, rather than (y+1)/U,, it is of interest to 
recalculate and replot the present results accordingly so as 
to ascertain to what degree these thoughts are borne out by 
actual experiment. The results are shown in figure 7. It 
can be seen that the theoretical and experimental results 
are in nearly perfect agreement. Comparison of figures 6 
and 7 provides striking evidence supporting the contention 
that k should be equated to M,(y+1)/U, rather than 
(y+1)/U,, as has been done so often in the past. 

Wings of finite aspect ratio.— The similarity rules for wings 
of finite aspect ratio are given by equations (59) and (60). 
Although no essential simplification of the rules occurs for 
wings of small aspect ratio, the range of applicability of 
linear theory increases as the aspect ratio decreases. ‘This 
point can be illustrated by considering the results provided 
by linear theory in a specific case. A good example to 
select for this purpose is that of the drag in supersonic flow 
of a triangular wing with symmetrical double-wedge airfoils. 
(See ref. 31.) This particular choice was made for the follow- 

321695—f5——67 
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Ficurs 7.—Drag correlation of single-wedge section, k= M3(y+1)/U,. 


-2 


ing reasons: (a) Solutions are known for all supersonic Mach 
numbers; (b) the double-wedge airfoil discussed in the pre- 
ceding sections corresponds to the limiting case of the wing 
of very great aspect ratio. The drag results provided for 
this wing by linear theory are presented in figure 8. 


oe T a/ f [Dy 


(t/c)? 


Me -| 


N- Puckett 
(Linear theory) 


I 
2 1e 
ME- Ana, A 


Fraure 8.—Supersonic drag of triangular wing, linear theory. 


The results of neue 8 are presented in a different manner 


in figure 9 wherein Op. is plotted as a function of & for 
various values of A as suggested by equation (59). For 
purposes of comparison, the curves for the drag of airfoils 
(= œ) computed by both linear and.nonlinear theory are 
also included in the graph. As noted in the preceding 
section on airfoils, comparison of the results of linear theory 
with those of nonlinear theory for wings of A= œ shows 
that good agreement exists for larger ¢, but that at smaller 


weer -— -7 — — - m ~ 


Nonlinear theory 


——— Linear theory 


o ] 2 3 
Eo 


Ficurb 9.—Variation with £, of reduced drag coefficient of triangular 


wing. 


t, the values predicted by linear theory become too large. 
For wings of finite A, only the results of linear theory are 
available. They also exhibit the trend of indicating infinite 
drag as $ apo zero; however, the range of & in 


which the value of A is excessive is much less than is the 


case for A=. In general, Ob, of wings of small aspect 
ratio diminishes with decreasing "aspect ratio. 
The drag results of figure 8 are presented in still oe 


form in figure 10 in which is plotted the variation of Ob. 
with A for various values of ~. The principal merit of 
this method of plotting is that it aids in distinguishing the 
region where nonlinear theory must be used from that where 
linear theory may be useful. Thus, the tsve-dimensional 
nonlinear theory results appear on the right of the graph 
corresponding to large A; whereas the three-dimensional 
linear theory results appear on the left for small A. The 
filling in of the remainder of the graph requires either the 
solution of the equations of three-dimensional nonlinear 
wing theory or the use of experimental data. It should 
again be noted that the present drag considerations apply 
only to the pressure drag. Before plotting experimental 


results in the manner indicated, it is necessary to first. 


subtract the friction drag. 
LIFT 
The similarity rule for the lift coefficient Cz of a family 
of wings having ordinates Z een by equation (44) 


Zle=t[e| le, y/B)— = | (44) 
is given in equation (46) as 
C= Ut) Cr=L (Eo, A, a) (46) 
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Nonlinear theory 


——-— Linear theory 


2 08) 
A 
FieurE 10.—Variation with A of reduced drag coefficient of triangular 
wings. 


where 


TOE 
[Uk 


If hypothesis (d) is accepted, Cz varies linearly. with «œ for 
at least small angles of attack. It is advantageous, therefore, 
to consider the lift ratio Cz/æ rather than Oz alone, thereby 
minimizing the influence of &. 

On __ pi (i/ ci Se, dl , 

a a (Uk)? LlEo, A, D= Toys L (Eo, 
where £’ (s, A, č) is anew function of the indicated variables 
obtained from £ (¢,, A, &) by division by &. Therefore, 
lift results may be presented by plotting the variation with 


A=[U.K(t/e)|"* A, iT, (49) 


A, &) (70) 


l £,, A, and & of a generalized lift ratio Ola defined as follows 


(the prime on £ has been omitted for simplicity): . 


N : 
C,Ja==[Uok(te)|"5 (Cr/a)=L (Eo, A, &) (71) 
= 

Equation (71) shows that C,/a depends upon three 
parameters, one more than the number which can readily be 
treated on a simple plot. Simplification can be gained, 
of course, by holding one of the parameters constant for an 
entire graph. Results so presented are particularly interest- 
ing for £=0, (M,=1); A= œ, (á= œ); and &=0, (a=0). 
The latter scheme is especially good since experiments 
indicate that lift curves of wings are often relatively straight 


——s 
lines at all Mach numbers. The values of C,/a at €=0 
might, therefore, be expected to be good indications of the 
actual values for other & The appropriate similarity rule 


"may then be written 


(Cale)=L to A, = Lato, A) (72) 
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For cases where linear theory applies (hypotheses (b) and 
(c)), two statements can be made immediately which 
provide further information about C,/a: (a) (C,/a): must 
be independent of k, and (b) (C,/a); must be independent of 
@ by virtue of the superposition principle of linear theory. 
Therefore, following the procedure used in equation (60) 
gives 


L: (maria), (2 = £ (VELA) 


r ië 
(73) 


where again £; is a different function for subsonic and super- 
sonic flow. 


(G-a 


Wings of infinite aspect ratio. or wings of infinite aspect - 
ratio, the functional relation of linear theory for the lift ratio 


given by equation (73) reduces to 


Te 


(2) _ const. _ 
g |M. JMe—i| 


(2) 


Solutions of the equations of linear theory show that the 
value of the constant is 27 for subsonic flow and four for 
supersonic flow. Examination of these results indicates that 
they are valid at Mach numbers appreciably less than or 
greater than unity, but are invalid for Mach numbers near 1. 

A. similarity rule for the section lift coefficients of a family 
of affinely related symmetrical airfoils which is valid through- 
out the Mach number range may be obtained from equation 
(71) by setting A= œ. 


(cy) =L (bo, œo , X= (fo 2a) 


At a Mach number of unity, is identically zero and the 
expression for the lift ratio becomes 


] a 
Cle) rE (0, a) (76) 


Equation (76), when considered together with hypothesis (d), 
indicates that at sonic speed, the lift-curve slope at zero 
angle of attack of airfoils of a single family varies inversely 
as the cube root of the thickness ratio, thus, 


(74) 


(75) 


S 
(c,/a) — Z(0 5a), 
E,=0 


const. 


[Uok (te) 


Note that as the thickness ratio goes to zero, the value of the 
lift-curve slope at zero angle of attack becomes infinite, just 
as is indicated by equation (74) to be the case according to 
linear theory. If, on the other hand, t/e diminishes to zero 
while æ is fixed so that @ becomes very large, it is plausible 
that the thickness ratio does not have any effect on Cr We 
thus have the following: 


ai 7 
(¢o/a) st =k or 1(0,a—-0)= 


6 4 ra ols x 
tel Uk (t/c)}*/ i (0, tfe —> 2 J= nX const. (78) 


In this case, c; is proportional to the two-thirds power of a. 


If hypothesis (e) is acceptable, the variation of Be with £, 
at M,=1 is zero, since 


(77) 
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eG A Em0 E ON! Em0 (9) 
-SOLO 


Tf k is taken to be Mè (y+1)/U., (dOpldMo)atem1 18 given by 
equation (56) and the variation of c,/a with Ms at M,=1 is 


ETORRI IRON, Paii 
; 3(=) | Af,=1 a(é t)=-3 ~ B8\a aN 


The corresponding exact relation can be determined similarly 
when equation (52) is used for the slope of the pressure curve. 


O hOla 


At present, calculations have been made of the lift in sonic 
flow (Guderley and Yoshihara, ref. 32) and in slightly super- 
sonic flow (Vincenti and Wagoner, ref. 10) of symmetrical 
double-wedge airfoils inclined an infinitesimal angle of attack. 
The corresponding solutions for airfoils in slightly subsonic 
flow have not yet been found. The results of the above- 
mentioned investigations are shown in figure 11 together with 


o a 81) 


Double— wedge airfoil y, 
Guderley and Yoshihara-—” ,/ 
Vincenti and Waqoner-——“ 
—3 =a -I O i 2 3 
3 


gJ 
Fraurs 11.—Theoretical lift-curve slope of double-wedge airfoils. 


the corresponding values given by linear theory. As in the 
drag case discussed in the preceding section, it may be seen 
that the remarks concerning the relation between the results 
of linear theory and nonlinear theory and the slope of the 


curve of c e versus & at &,=0 are verified for this particular 
airfoil. 

Wings of vanishing aspect ratio.—T wo well-known results 
of linear theory are that the lift-curve slopes of wings of 
finite aspect ratio remain finite throughout the entire Mach 
number range and that the lift-curve slopes of wings of 


vanishing aspect ratio are mdependent of Mach number, 


— 

Therefore, equation (73) implies that (C,/a); must be pro- 
portional to 4 either for wings of vanishing A in any flow, 
or for any wing in a flow of vanishing £: 


(Cr/e)i=(CrJa=A Xconst., (Cz/a)i=(Cz/o)1=-A X const. (82) 
A—/) fad Ad Mex} 


The value of the constant must be determined for each plan 
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form by actually solvmg the equations of linear theory. 
For wings having trailing edges which possess no cutouts 
extending forward of the most forward station of maximum 
span, that is, triangular, rectangular, elliptical, etc., as well 
as certain sweptback wings, the value of the constant is 
aw/2. (Refs. 26 through 29 should be consulted for further 
discussion of this point as well as for the values of the con- 
stant for wings having cutouts m the trailing edge which 
violate the above-stated condition.) 

It is seen from equation (82) that the lift-curve slopes of 
wings in sonic flow decrease continuously in magnitude as 
the aspect ratio diminishes toward zero. Since, in addition, 
the lift-curve slope given by linear theory has its maximum 
value at M,=1, it is conjectured that the lift results of linear 
theory are a good approximation to those of nonlinear 
theory not only for all wings at Mach numbers far from unity 
but also for all Mach numbers for wings of sufficiently small 
aspect ratio. 

Wings of finite aspect ratio— At the present (iae no solu- 
tions of the nonlinear theory are available for wings of finite 
aspect ratio. However, from the remarks of the preceding 
paragraphs, it is apparent that a curve representing the varia- 


N as 
tion of C,/a with A for constant & and & would have the 


following asymptotic properties: Cr would increase linearly 
with A for small A and be independent of A for large A. 
In order to give a better idea of the numerical values to. be 
expected, a set of typical results of this type is shown in 
figure 12 for wings having triangular plan forms and sym- 


J Nonlinear theory, @ =O 


— —— Linear theory 


O 2 4 
A 
Figure 12.—Variation with A of reduced lift-curve slope of triangular 
i wings. 


metrical double-wedge airfoil sections. The supersonic 
results are those of Stewart, Brown (refs. 33 and 34), and 
others. The subsonic results are those calculated by De 
Young and Harper (ref. 35) using Weissinger’s modified 
lifting-line theory. 


ee 
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An interesting result of the foregoing remarks concerns the 
influence of thickness ratio on the lift-curve slope at zero 
angle of attack of wings in sonic flow. For wings of large 
aspect ratio, the lift-curve slope is inversely proportional to 
the cube root of the thickness ratio. . For wings of small 
aspect ratio, the lift-curve slope is independent of the thick- 
ness ratio. 

PITCHING MOMENT 


The remarks of the pitching-moment characteristics of 
wings follow in a manner exactly analogous to those just 
stated for the lift characteristics. The corresponding state- 
ments for the pitching-moment coefficient C,, may be obtained 
by simply replacing Cz with C, and £ with M. Thus, the 
similarity rules for C,, corresponding to equations (71) and 
(73) are the following, respectively: 


(Ol = [UEN Cn) = (to, A, 2 

(ee) ay A), 
Cr =l MAET) (84 
( =) Te] (Viil 4) (84) 


(83) 


where once more M and M, are different functions for sub- 


sonic and supersonic flow. The only difference between the 
discussion of Cu and C; is that the values of the constants of 
equations (74), (77), (78), and (82) are, of course, different. 
Graphs of theoretical pitching-moment characteristics for 
airfoils and for triangular wings corresponding to the lift 
results of figures 11 and 12 are shown in figures 13 and 14. 


Double—wedge airfoil P 
Guderley and Yoshihara—~ 
Vincenti and Wagoner |——-” 


Nonlinear theory, a» 0 
~——— Linear theory 


Figure 13.—Theoretical pitching-moment curve slope of double-wedge 
airfoils. 


The moment axis is taken to be through the most forward 
point of the wing. 

Sometimes it is desired to present pitching-moment char- 
acteristics of wings in terms of center-of-pressure position 
rather than pitching-moment coefficient. Since the center- 
of-presspre position can be expressed in terms of Om and C, by 


Tep. Cn 


— 5 


r O, (85) 
the resulting expression for the center-of-pressure position 
found through application. of equations (71) and (83) is 


on Min AD CP En ED 


86 
© Kén AD i 
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Nonlinear theory, a = O 


—— — — Linear theory 


0 g 4 œ% 


M 


A 
Figure 14.—Variation with A of reduced pitching-moment curve slope 


of triangular wings. 
The corresponding relation for linear theory is 
(4222) =C. Pa WIA) =C. Pal) (87) 


PRESSURE DRAG DUE TO LIFT 


The similarity rule for the pressure drag of inclined wings ` 


having symmetrical airfoils is indicated by equation (48) to 
be the following: 


Ok ~ 
o= eiar Co= Dln A, 2) (88) 
The portion of the drag due to lift is therefore 
ACp=Co—Co,=crp ER ae [D(t,, A, 2)—D¢,, A, 0)] 
= fia Dalin Z A, @) i (89) 


Since the differential pressures between upper and lower sur- 
faces of wings having symmetrical airfoil sections are pro- 
portional to æ for at least a small range of œ surrounding zero 
(bypothesis (d)), it follows that ACp is proportional to the 
square of a. It is advantageous, therefore, to consider the 
drag-rise ratio ACp/a? rather than AC> alone, thus, 

OI YO Dilly LO rE PEA) (00) 
where D,’ is a new function of the indicated variables ob- 
tained from Da by dividing by a. Consequently, drag-due- 
to-lift data should be presented by plotting the variation with 


Eo A, and @ of a generalized drag-rise ratio ACp/a* defined as 
follows (the prime on “Dy, is omitted for simplicity): 


AC plat=[Upk(t/c)]/{ACp/o2)=Dalko A, &) (91) 
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The actual presentation of the results of this three-parameter 
system may be accomplished as described in the section on 
the lift of wings. ° Of particular interest is the simplification 
resulting from presenting only the values found at a=0. 
The simplified similarity rule is then 


CO nfo) =Dalbn A, 0) (92) 


For cases where linear theory applies, the following results 
hold: 


(AS) - Ta il (FTI A), 


CH) 


Wings of infinite aspect ratio.— For wings of infinite aspect 
ratio, the functional relation of linear theory for the drag due 
to lift, equation (93), reduces to 

const. Se) 


aa? 


Solutions of the equations of linear tice show that the 
value of the constant is zero for subsonic flow and four for 
supersonic flow about any symmetrical airfoil. These results 
are valid at Mach numbers appreciably less than or greater 
than unity but are invalid for Mach numbers near 1. 

A similarity rule for the drag due to lift of a family of 
affinely related symmetrical airfoils that is valid throughout 
the Mach number range may be obtained from equation (91) 
by setting A= œ. 


Fey Po (Me A) (93) 


= 


(94) 


(ŽE) = Dal, o, B= dal bo, 2 (95) 


At a Mach number of unity, equation (95), for the drag due 
to lift, reduces to the following: 


ae l 
OE” aa( 0 i) (96) 


Equation (96), together with hypothesis (d), indicates that, 
at sonic speed, the drag-rise ratio Acg/a* of airfoils of a single 
family varies inversely as the cube root of the thickness ratio. 
For very large values of @, the thickness ratio cannot have 
any effect on Aca; therefore, Aca is proportional to the five- 
thirds power of the angle of attack. If hypothesis (e) is 


a 
accepted, the variation of Acg/a’ with & at 17,=1 is zero, 


CM, ZL @..]¢@ 
-0 & a (Ea? E) 


If k is taken to be M? (y+1)/Uo, (€dC,/dM_)arj1 is, given by 
equation (56), and the variation of Acg/a? with M, at 14,=1 is 


C =d, =: 


(97) 
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EAA Pade | ate bn’ G) 
=). Fal ar (2 Jaa? G -a a we 


The corresponding exact relation can be determined similarly 
by use of equation (52) for the slope of the pressure curve. 


Itis 
[ar (s aa r a 


Af o—> 


Wings of vanishing aspect ratió.—Since, as`indicated in 
the preceding sections, the lift calculated by means of lmear 
theory remains finite throughout the Mach number range, 
and diminishes to zero as the aspect ratio approaches zero, 
it follows that the drag due to lift must behave similarly. 
It is consequently reasonable to presume that linear theory is 
capable of describing the drag-due-to-lift characteristics of 


wings of vanishing aspect -ratio at all Mach numbers. 


Therefore, the following relations stemming from equation. 


(93) hold: 


CA, -6 ») =A Xconst., 


(29) - B 


ac D 
Solutions of the equations of linear theory show that the 
value of the constant is 7/4 for all wings of small 4/|A7,2—1].4 


(100) 


= A const. 
=] 


n tee o o o a am = 
+ 
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whose trailing edges possess no cutouts extending forward of 
the most forward station of maximum span (i. e., triangular, 
rectangular, elliptical wings, etc.). The quoted value of the 
constant corresponds to the development of the full “leading- 
edge force.” It is known -that this force is oftentimes not 
completely realized, due to a local separation and subsequent 
reattachment of the flow around the leading edge. If the 
leading-edge force is nonexistent, the corresponding value 
for the constant is 7/2. 

Wings of finite aspect ratio —At the present time no drag- 
due-to-lift results have been obtained from the nonlinear 
theory for wings of either finite or infinite aspect ratio. 
The foregoing remarks, however, are sufficient to determine 


that a curve representing the variation of ACp/a* with A 


for constant £, and & would increase linearly with A for 


small A (unless the degree of attainment of the leading-edge 
force also depends on A) and become independent of A for 
large A. The resulting curve would presumably have the 


same general appearance as that shown in figure 12 for Ch Joe 

It may sometimes be desired to present drag-due-to-lift 
results in terms of AC>/C;,? or AC>/aC, rather than ACp/c’. 
The similarity rules for these quantities can be quickly de- 
duced from the foregoing results. 


Ames ABRONAUTICAL LABORATORY 
NATIONAL ÅDVISORY CoamuTTEeE FOR AERONAUTICS 
_ Morrerr FreD, Carr., Dec. 16, 1962 


APPENDIX A 
RELATION BETWEEN U, AND a* STATEMENTS OF THE TRANSONIC-FLOW EQUATIONS 


Equations (3), (7), (12), (13), (15), and (17) were presented 
in the text as being applicable to the study of transonic, as 
well as supersonic, flow about thin wings. These equations 
repeated below as equations (Al) through (A6), will be 
referred to as the U, statement of the problem since the 
perturbation velocities are taken about the free-stream veloc- 
ity U. The perturbation potential ¢ is defined by 


g=—U 2+ (A1) 
The differential equation is 
(1 — M,?) Pret Pyy t Crh Prr (A2) 
The shock relation is 
(1 ~ MM) (2z — 2z) + (2r — Pn) + (on — Pa) 
Pa Fe 
=k (7) een) A3 
The boundary conditions are 
ab r=— o 
a (yz) o— (9y) gm (p:) o=0 (A4) 
at the wing surface 
j (P:)z=0= Uo (£) (A5) 
The pressure coefficient is given by 
2 
Cp <= FT U, Px (A6) 


In the usual a* statement of the transonic-flow equations 
(eqs. (24) through (28) in the text), it is assumed that all 
velocities are only slightly different from the critical speed 


of sound a*. The perturbation potential is defined by 
gp =—a*z+h (A7) 
The differential equation is 
; l , 
Pwt p'am 0" x00" ne (A8) 


The shock relation is 


(g ye y +(¢’ ry 9 y= 


yrl (£ nt? x CAR (AQ) 


a” 2 


If the perturbation analysis is carried out in a completely 
consistent manner, the boundary conditions are: 
at r—=— © 


aC ‘o= 


aie 1’ (o’ o=(¢'s)o= 0 (A10) 


at the wing surface 
z ) — 7* OZ A 1 
(Y z z=0™ 4 Ox ( 1 ) 
and the pressure coefficient is given by 
ee 
p= — lo'le (A12) 


The relation between the U, and a* statements can be 
determined directly in the following manner. The differ- 
ential equations and shock relations for ¢’ and ¢ will be the 
same if 


1 òg’ 
rh De Peay 
or if 
sr p © —_a* (1— M5) 
oak y+1° y+1 ey 


The boundary conditions for »’ corresponding to those 
stated for y in equations (A4) and (A5) are: i 


at t= — © 
EOE N _a* (1—, 0) _a* (1—M,*) (1—M,") 
(g zo=k y-+1 (CAF yI Jel 
a* a* (A14) 
(a=k Tz C50, (pk Tz (¥s)o=0 
at the wing surface 
Iy k a* k a* A 5 
(y dai +1 (¢2)r20= Uo y+1 (£ ( 15) 
Finally, the pressure coefficient is given by 
? — 2 Po fof as 2k kU, 
ee Se ae HO 


Comparison of the above equations with those given 
previously for the completely consistent a* analysis reveals 
their identity if, and only if, k is taken as (y+1)/U,, that is, 


_ytl 
Uo 


(A17 


As far as obtaining the values of C, is concerned, therefore, 
the conventional a* statement of the problem may be 
regarded as a transformation of the more general U, state- 
ment with & being defined as in equation (A17) or (19). As 
is evident from equation (A13), however, the local velocities 
found in the a* analysis are only correct when the free- 
stream Mach number is unity. 
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The foregoing discussion has demonstrated that the con- 
ventional a* statement of transonic-flow problems yields 
the same values for C, as the U, statement, provided k in 
the latter set of equations is equated to (y+1)/U.. It is 
indicated in the text, however, that transonic results ob- 
tained by use of other expressions for & (and, in particular, 
k= (y+1)/U.) are superior in many important aspects 
to those obtained using k=(y+1)/U.. It is consequently 
not without interest to determine a generalized form of the 
a* equations which correspond to the U, equations with 
unspecified k. The resulting equations for the perturbation 
potential and the shock relation are: 


(A18) 
(A19) 


p'=— a*r 49 
Pyt p a=k p 10's 
/ t 
(0! — Ow) HE ae a Ek ois aii =) (o'a P a (A20) 


The boundary conditions are: 
at r= — œ 


Mè— 1 
po’ 


(y= C ’ Do =0 (A21 ) 


(g'z) = 


REPORT 1153— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


at the wing surface 


sa a* 65 


(A22) 
and the pressure coefficient is given by 
, 2 o , 

Cr =i [p's —(¢’s)ol (A23) 


Tt can be readily verified that these equations will produce 
the same values for the pressure coefficient as equations (A1) 
through (A6), provided k and k’ are related as follows: 


T aF (A24) 


Thus, if it is desired to obtain values for O, by use of equa- 
tions (A18) through (A23) that are the same as those given 
by equations (Al) through (A6) with k=Mè(yt+1)/U., k’ 
should be equated to the following: `“ 


pr Mo (y+ 1) 


a 


(A25) 


APPENDIX B 


DERIVATION OF SIMILARITY RULES 


The basic equations of linear N and of nonlinear theory 
of transonic flow may be summarized as follows. The 
differential equations are: 


(B1) 
(B2) 


42 Si, 0 Linear 

(1— -MA le £ Sar eae mee 

a oy ae k— Be Og Nonlinear 
oz Dr? 


The approximate shock relation used in the treatment of 
transonic flow is 


(1—M,") (91 — a) + (4, — 997)? + (enen) 


=1( 21%) (Pa — Oxy)” (B3) 
The boundary conditions are: 
at r=— o 
Dp - 
Ionon. wa 
at the wing surface 
Le) 3,9 _ 7(2 4 
U, 9 ajo) I (= ; (Re 
where the geometry of the wing is given by’ 
(Z/c)=7f(@/c, y/6) (B6) 
The pressure coefficient is given by 
20 
a= =U, T (B7) 


If the differential equations are now transformed into a 
system with primed quantities and the proportionality or 
stretching factors are denoted by s with appropriate sub- 
scripts such that 


U,/=sy Uo 


= S820, y= Sy; 2! = 8,2, y= 89, 
aK m 8 (B8) 
y (1 — M.) =88 1— M. 3 Uk’ =8, Uk or kik 
U 
the potential equations (eqs. (B1) and (B2)) become 
e Be e Die ot ie 
Or?” 35 OYT 83:02" 
0 Linear (B9) 


Sule ,, 09’ 070’ ; 
EE k TETE Nonlinear (B10) 


and the shock relation becomes 


321605—55 68 


3.2 
A (1 Paia (ovi 2 ya) *+ 


= sete k’ Lan 


E Ae z’ —p': 3) aes (o'z e e) 


Nonlinear (B11) 


The flow in the primed system is similar to that in the 
original system if y’ satisfies the same differential equations 
and boundary conditions as y. Consequently, for similarity 
to exist, it is necessary first of all that the potential equations 
and the shock equations for the two flows be the same. 
Therefore, the following relations between the stretching 
factors must be satisfied: ose 

ZN 8v j £ os Linear. (B12) 

Sz Sz 8x80 (s,7/s, Nonlinear (B13) 


where, for linear theory, A/N is an arbitrary constant which 
can be equated to 8,*/s, if desired. The constant is written 


‘as a fraction in order to maintain a certain symmetry 


throughout the analysis. 

An immediate consequence of this transformation is that 
the wing plan forms undergo an affine transformation such 
that the aspect ratios of wings in similar flow fields are 
related, according to both linear and nonlinear theory, by 


,_8 4_1 1—M.? 

Aint A= A= aay 4 (B14) 
or by 

V-M A’=J1—M? A (B15) 


Since g’ is proportional to y, the boundary conditions at 
g=—o are automatically satisfied. The boundary condi- 
tions at the m may be given in either of two forms: 


og" &— (0 __ Se Ò (2 x) 

€ ee oz yao. -8e Uor JODA c b (B16) 

Oy’ To Ò ; z y’ m , ò (z y' 

D2 armo T IET (7 yy) 800. EEDA P) 
(B17) 


whence, if the two wings have the same ordinate-distribu- 
tion functions, that is, if f (æ'/ce, y'[b)=f(ele, y/b), the 
ordinate-amplitude parameters are related as follows: 


1—8? $088 

Su8: Su 8z 
à A Linear (B18) 
Er Ger | (aie T Nonlinear (B18) 
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(OMY Mn 
a (ae 
vV — Me) — M7)’ _y1—M 


SNE 
UA TIG RUADA Nonlinear (B21) 


The relation between the pressure coefficients at corre- 
sponding points on the wing surface is given by 


Olea. Se) _ 2 ðe) _ 8 
ý N 1'0 ere UJ, 0x) imo Sus: ” 
(x Ca Linear (B22) 


fr 23/7 Uk 1/3 . 
=) (GE Cy Nonlinear (B23) 


or more completely by 
(ary y! 
Cy e 7 oy, -MAF o) A’; = zg] 


“po (ER mai Linéar AB) 


Cf Rovers JMF 4; 5 n% 


(7) CY @ Te 
Ji—M A; zg] Nonlinear (B25) 


The foregoing relationships may be summarized in ‘the 
following statement: The similarity rule for the pressure 
coefficients on a family of wings having their geometry given 


by 
(Z/c)=7f(@/e, y/b) (B26) 


1S 
= 2- 
=l P, (a Vi-M} A; =, 9 Linear (B27) 


C= GR: ae? | es (jkr = aie FASS Nonlinear (B28) 


The similarity rules for the lift, pied E and drag 
coefficients given by the linear theory are : 


> 


af za 2 : 
Cn => Ly oo V1i—M? 4) Linear (B29) 
Cat M: (vine Vin A). Linear (B30) 


cme 2 
Cry=2 D, as VI-M; 4) Linear (B31) 


The corresponding similarity rules given by nonlinear theory 
are 


2/3 a We 
C= ryn L Eaa q I — M a| Nonlinear (B32) 


cn ce te | 


(B20). 
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C= TI n M pi (Ukr y1— MF A | Nonlinear (B33) 


Co= ‘ai ore D oles Vay VI— TT A | Nonlinear (B34) 


It should be noted that the foregoing equations have been 
written for subsonic flow where M,<1. If M,>1, the 
radical -/1—M,? should be replaced with ¥M?—1. 

In the linear-theory analysis, ` has remained a completely 
arbitrary coefficient to be selected as best suits the particular 
problem at hand. For instance, the compressible-flow rela- 
tionships between two wings having identical pressure dis- 
tributions are found by setting A=1. If, on the other hand, 
à is set equal to y1—M ê, the thickness ratio, camber, and 
angle of attack of related wings are identical. The greatest 
simplification of the similarity rules of linear theory occurs 


when à is set equal to -(1—M/,2/r since the number of param- 
eters necessary to show the results of linear theory is thereby 
decreased by one. This degree of arbitrariness in the 
similarity rules for linear theory is in contrast to the case for 
nonlinear transonic theory in which no undetermined co- 
efficient like à is to be found. 

For the present purpose of gaining a better understanding 
of transonic flows, the most advantageous choice for A is 


oft 1/3 
eee (B36) 


because then the similarity rules for linear theory assume 
forms identical to those for nonlinear transonic theory. This 
is important since it implies that solutions of linear theory 
and of nonlinear transonic theory can be expressed’ as func- 
tions of the same parameters and, hence, can both be plotted 
on a single graph in terms of one set of parameters. The two 
theories would, of course, yield two distinct curves on such a 
plot. The curve for linear theory would be accepted as valid 
for purely subsonic and purely supersonic flows but may or 
mey not be valid in the transonic range, as discussed pro- 
viously. The curve for nonlinear transonic theory is valid 
not only for transonic flows, but for subsonic and supersonic 
small perturbation flows as well. As is evident from the 
derivation of the basic equations, however, the results of the 
nonlineartransonic theory should be considered to be of only 
equal accuracy to those of linear theory in the definitely 
subsonic and supersonic regimes, despite the fact that the 


solutions are much more difficult to obtain. 
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